July 1969 


Brief 69-10204 



AEC-NASA TECH BRIEF 




AEC-NASA Tech Briefs describe innovations resulting from the research and development program 
of the U.S. AEC or from AEC-NASA interagency efforts. They are issued to encourage commercial 
application. Tech Briefs are published by NASA and may be purchased, at 15 cents each, from 
the Clearinghouse for Federal Scientific and Technical Information, Springfield, Virginia 22151. 


Some Numerical Methods 
of First-Order Ordinary 

Study of some numerical methods of integration of 
systems of coupled first-order ordinary differential 
equations has been reported (/). The report describes 
the extrapolation methods of Bulirsch-Stoer and 
Neville, presents the results of tests comparing these 
extrapolation methods with the Runge-Kutta and 
Adams-Moulton methods, and shows some circum- 
stances under which the extrapolation method may be 
preferred. Results of double-precision tests of the 
Bulirsch-Stoer method show the efficiency of various 
orders of this method as a function of desired accuracy. 

The speed and accuracy of the methods were com- 
pared by programming the methods on the same 
computer (a CDC-3600) and applying each to the 
same sets of equations. The numbers of times that 
the functions to be integrated had to be evaluated, for 
a desired accuracy, were used as a computer-indepen- 
dent and program-independent measure of time. The 
accuracies were compared by use of each method to 
integrate from given initial conditions to some final 
point at which the relative error was computed. 

The results of cases run by the four methods tested 
indicate that the efficiency of the solution is quite 
dependent on proper matching of the order of the 
method of solution to the desired accuracy. For ex- 
ample, in tests on the negative-exponental equation, 
for a relative error of 10 - lf) the Neville method with 
M =4 takes twice as long as the method with M =6 and 
about 2.5 times as long as the shortest time with M = 10. 
The Runge-Kutta fourth-order method and the Adams- 
Moulton fourth-order method take about 160 times 
as long as the Neville method with M - 10, or about 
250 times as long as the Bulirsch-Stoer method with 
M - 10. Even at relative errors of 10“ J the Adams- 
Moulton and Runge-Kutta methods take about 9 times 
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as long as the Neville and Bulirsch-Stoer methods 
with M =6. At the other extreme, if a rough estimate 
is desired and errors up to 10“ 2 can be tolerated, 
Neville's method with M = 10 takes almost 7 times as 
long as the fourth-order Adams-Moulton methods, and 
the methods of lower order are preferable. 

The double-precision tests of the Bulirsch-Stoer 
method show it to be quite powerful for work of great 
accuracy; it seemed especially capable of coping with 
changing difficulty bf solution of the equations pre- 
sented to it. 
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